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The paper is devoted to vector or pseudoscalar heavy quarkonia production ip( \ Bs in as- 
sociation with five charged 7r-mesons in exclusive _B c -meson decays. Using available formfactors 
parameterizations and spectral functions of W — > 5n transition we obtain branching fractions of 
these decays and distributions over invariant mass of (57r)-system. 
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I. INTRODUCTION 



o 

^jj^ Heavy quarkonia (for example charmonium mesons rj Cl J/?/ 1 , ip(2S), etc, or bottomonia rjb, T, Xbj, ■ ■■) always 
■ played a special role in elementary particle physics. Due to the presence of a heavy quark with mass toq 3> Aqcd the 
scales of quark-antiquark annihilation (~ I/tuq) and hadronization into experimentally observed meson (~ 1/Aqcd) 
differ significantly and these two processes are separated. As a result the reactions of heavy quarkonia production 
and decays can be used for analysis of strong interaction both in perturbative and nonperturbative regimes. Over the 
last years significant theoretical and experimental results were obtained in this field. 

Heavy quarkonia with open flavor (6c), i.e. _B c -meson and its excitations, take intermediate place between char- 
Cl^ monium and bottomonium mesons. So, they provide a possibility for independent test of models used in analysis of 
JL. . states with hidden flavor. Theoretical predictions for the width and lifetime of B c -meson[TJ 

at! 

M Bc = 6.25 GeV, r Bc =0.46ps (1) 

are in good agreement with experimental values [2|. Predictions for partial widths, on the other hand, differ signifi- 
cantly from experimental results. For example, ratios presented in ref. Q 

00 

OS : cr B Br(B c ^J/^) = 

CO ■ creBr (B -> J/ipK) y ' 



and 



——^——j—^- 0.249 ± 0.045_ . 076 (3) 

are approximately an order of magnitude higher than theoretical expectations based on current estimates for £? c -meson 
, production cross section and the branching fraction of its semileptonic decay 

In the present paper we consider exclusive decays B c — > V(P) + TZ, where TZ is a set of light mesons (e.g. 5n). 
According to factorization theorem the widths of these decays are connected directly with the widths of r-lepton 
decays r — > v T + TZ. In both cases the final state TZ is controlled by virtual VF-boson that is produced in heavy quark 
weak decays for hadronic processes or r — > v T W in leptonic ones. In contrast to r-lepton decays, B c — > V(P)W 
transition is described by formfactors, so it is possible to determine these formfactors using available experimental 
data about r — > v T TZ and B c — > V(P)TZ decays and compare them with theoretical predictions based on QCD sum 
rules and different potential models. Analysis of B c — > V(P)TZ decays can also be used to describe W — > TZ transitions 
in energy regions that cannot be achieved in r-lepton decays. 

In the next section of our paper we present analytical expressions for B c — > V(P)TZ widths and distributions over 
invariant mass of light meson system TZ. Technique of spectral functions used in our paper is also described briefly 
in that section. In section [TTT1 we give numerical expressions for B c — > V(P) fromfactors, W — > TZ spectral functions 
and present our predictions for branching fractions of decays and q 2 distributions considered in our article. Brief 
discussion of the results is given in the last section. 
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Figure 1: Typical diagram for B c — > V(P) + nn process 



II. ANALYTICAL RESULTS 



It is well known that in the valence approximation £? c -mesons consist of b- and c-quarks, so their decays into heavy 
quarkonia are caused by weak decays of one of the constituent quarks: b — > cW* for charmonium meson in the final 
state or c — > sW* in the case of Bg production. Produced in this decay heavy quark in combination with other B c 
constituent forms final heavy quarkonium, while virtual M^-boson hadronizes into light mesons system 1Z. Typical 
diagram of this process is shown in figfT] 

In our paper we consider decays 



B c -> V(P)+TZ, (4) 

where V(P) is vector (pseudoscalar) heavy quarkonium J/tp, tp(2S) or Bs , and 1Z = (57r) c his a set of five charged tt- 
mesons (production of smaller number of 7r's was considered in previous articles The amplitude of the reaction 

(j4|) in the framework of factorization theorem can be expressed through formfactors of B c — > V(P) transition and is 
equal to 

M V{P) = ^i^eW, (5) 

where Vu is the element of CKM mixing matrix, e^ w > is the effective polarization vector of virtual H^-boson, coefficient 
cl\ describes the effect of soft gluon rescattering jjj, and Hycp) is the amplitude of B c — > V(P) transition. This 
amplitude can be written in the form 



^(<? £ ) , i Tit \ a („2\ ( jh ^il e ) 



H v = 2M 2 A (q 2 ) ^fi + ( Ml + M a ) A 1 (q 2 ) 



a am (ge) f/ i / M ?- M L A 2 ^(g 2 ) c Re ^B () 



Ah + M 2 V q 2 J Mi + M 

for vector meson production and 



< = n OT (,f + A - + Fo V) Mi-Mi,, (7) 

in the case of pseudoscalar meson in the final state. In the above expressions pi.2 and M\p are momenta and masses 
of initial and final heavy quarkonia, e M is the polarization vector of final vector meson, q — p± — p 2 is the momentum 
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of virtual VF-boson, i*o,+ (q 2 ), -^0,1.2 (q 2 ) an d V (q 2 ) are formfactors. It is clear that these formfactors cannot be 
determined from perturbative QCD, so one should apply some nonperturbative methods. In our paper we use the 
results based on QCD sum rules Q and solution of different potential models 0-EIJ, [ll||. Further these formfactor 
sets will be denoted as SR, PM1, and PM2 respectively. 

If we are interested only in partial widths of the considered decays and g 2 -distributions, it is convenient to integrate 
over the phase space of light mesons system and use the technique of spectral functions (see ref.[l2j for detailed 
description). In the framework of this approach the differential widths of _B c -meson decays into vector or pseudoscalar 
quarkonium have the form 



and 



dT [B c -> VIZ] 



G 2 F V 2 h a 2 M\\V n 



128ttM| (Mi + M 2 ) 



r/? 3 X 



(Mi + M 2 



12M 2 V 



M?f3 2 



A\ + Mfj3 2 Aj + SMgq^V 



2„2xr2 



2 (Mi + M 2 ) 2 (Ml - Mf - q 2 ) AiA a ] + AA 2 (Mi + M 2 f p?} , 



dT[B c ^PK] _ G 2 F a 2 \V i3 Y 



■p{(Ml - M 2 ) 2 F 2 P f + MtP 2 F + p*} . 



dq 32ttMi 

respectively. In these expressions j3 is the velocity of final quarkonium in _B c -meson rest frame 



(8) 



(9) 



= 



\ 



1 - 



M 2 + y/ q 
Mi 



\ 



1 - 



Mg - y /g : 
Mi 



and p£t \1 ) are longitudinal and transverse spectral functions defined according to 



(10) 



{q^u-q 2 9^) P T{q 2 )+q^PL{q 2 ) = ^ j d$ n (q -+ K) S ^ [e^f , 



where 



^ / (2ir 

is the Lorentz-invariant phase space of system 1Z. 



d?h 

3T 



(2^) J IE, 



(11) 



(12) 



III. NUMERICAL RESULTS 

Let us first consider explicit expressions for spectral functions of different final states that will be used in our article. 
The amplitude of W — > tt transition has the form 

^(k^J^W) = Uk», (13) 
where f„ ~ 130 MeV. Resulting transverse and longitudinal spectral functions are equal to 

(q 2 ) = 0, (q 2 ) = f 2 . (14) 

If we have five charged 7r-mesons in the final statefi.e. TZ = (57r) ch = 7r + 7r + 7r + 7r~7r~) the amplitude of W — > (57r) c h 
transition can be obtained from resonance model [13| (see diagram shown in fig'H]). Due to the partial conservation 
of vector current the longitudinal spectral function is equal to zero, while the transverse one can be approximated by 
the expression 
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(a) 



32 



25m" 



1 - 1.65s + 0.69s 2 



(s + 2.21) -4.69 



(15) 



where s is measured in GeV 2 (the spectral function itself is dimensionless). The dependence of this function on q 2 is 
shown in figEJ 

It should be noted that in addition to direct decays J/tp and five charged 7r-mesons can be produced in the reaction 
B c —> ip(2S) + (37r) c h —> J/ip + (57r) c h, so the spectral function of (37r) c h state should also be considered. As in 
the previous case, this transition can be described in the framework of resonance model with diagrams shown in 
figGIa) ( see ref.[lj|)- The contribution of longitudinal spectral function can be neglected, and the transverse one is 
approximated by the expression (see figEJb) ) 



s - 9m 2 \ 4 1 + 190s 

s ' [(s- 1.06) 2 + 0.48j 2 

Experimentally these two processes can be easily separated. 

Formfactors of i? c -meson transition into vector or pseudoscalar quarkonium can be calculated, for example, in the 
framework of QCD sum rules or various potential models. In our work we use formfactor sets presented in papers 
(denoted hereafter as SR), [j-fioj (PM1) and Hi] (PM2). The values of these formfactors at points q 2 = and 
q 2 = <7^ ax = (M 2 — M 2 ) I (2M\) are given in tables HI and ITT1 It can be easily seen from these tables that for ground 
quarkonium states formfactors' behaviour for different models is similar, while predictions for excited ?/)(2<S')-meson 
of PM2 model differ dramatically from others. For example, PM2 model formfactors decrease with the increase of 
the squared transferred momentum, while formfactors from SR and PM1 models increase. Such difference can be 
probably explained in the following way: in the framework of potential models _B c -meson formfactors are determined 
from the overlap of initial and final quarkonia wave functions; the wave function of ■0(2S')-meson (in contrast to ground 
states J '/if} and B^) has a node that leads to such unusual behaviour of the formfactors. It will be shown later that 
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Table I: Formfactors of B c — > V transition 
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Table II: Formfactors of B c — > B s transition 



such effect leads to large difference between different models' predictions for (^-distributions in B c — > i/j(2S) + (57r) c h 
decays. 

Numerical values of a\ coefficient for B c — > ipO + TZ and B c — > + 7?. decays are 



and 



a\ (m c ) = 1.14 
ai (mfc) = 1-2 



(16) 
(17) 



respectively. 

Substituting these numbers into relations ([5]) and (0) it is easy to obtain numerical values of the branching fractions 
of the decays considered in our article (see table Uri]) . In order to compare them with the experimental results it is 
also useful to consider the ratio of B c — > V(P) + 5ir and B c — > V(P) + ir branching fractions, where the dependence 
on the choice of a formfactor model is partially canceled. In the case of B c — > J/ip + TZ decay these ratios for different 
models are equal to 



Br [B c J/tp + (3tt) 



chJ 



Br [B c ->• J/ipn] 
that agrees well with experimental result [15| 

Bv exp [B c J/ijj + (3n) ch \ 



= (2.3) SR , (2.2) pM1 , (2.1) 



PM2 ' 



Br exp [B c -> J/tpn] 



= 2.4±0.3±0.3. 



and 



For other decays we get 
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(18) 

(19) 

(20) 
(21) 
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Table III: Branching fractions of B c — > V(P) + (57r) c h decays 




Figure 4: Tranferred momentum distributions for B c — > J/tjj + (37r) c hand B c — » ip(2S) + (37r) c h decays (left and right figures 
respectively). Solid, dashed and dotted lines correspond to formfactor sets SR, PM1, and PM2 



where in the last case we take into account contributions from both vector and pseudoscalar £? s -mesons. One can 
notice some interesting properties of the ratios presented above. First of all, branching fractions of charmonium 
production in association with one and five 7r-mesons are comparable with each other, while for £? s -meson in the final 
state the production of large number of 7r-mesons is strongly suppressed. This effect was observed also for 4-7T final 
state (see ref.[H) and is caused by the small difference between B c - and Bj-'-meson masses. In addition, one can 
see that PM2 prediction for B c — > ijj(2S)7Z decay is about an order of magnitude smaller, than SR and PM1 results. 
This difference can be explained by the behaviour of PM2 model formfactors for excited charmonium state mentioned 
above. 

Using relations © one can also obtain (^-distributions for the branching fractions considered in our article. 
In the case of B c — > J/ip + (37r) c h and B c — > ip(2S) + (37r) c h decays these distributions are shown in figHJ Available 
experimental results from [15J are shown in the left plot of this figure with dots. One can see that theoretical 
predictions based on SR formfactor model are in good agreement with the experimental data. 

Similar distributions for B c — > i/jO + (57r) c h and B c B^ + (nir)^ decays are shown in figsfS] and Irrespectively. 
One can easily see that for B c — > J/tp + (57r) c h the results of different formfactor models agree with each other up to 
overall normalization. In the case of B c — > i/j(2S) + (57r) c h decay, on the other hand, PM2 model prediction differs 
strongly from SR and PM1 results. This difference is caused by the behaviour of the formfactors mentioned above. 



IV. CONCLUSION 



B c -meson, i.e. heavy quarkonium build from b- and c-quarks, was discovered by ALEPH collaboration in 1997 [16| . 
Reported values of its mass and lifetime are in excellent agreement with theoretical predictions [l| . The situation with 
branching fractions of different decays (e.g B c — > J/ip£v or B c — > J/ipir) is much worth. It is clear that additional 
investigation of this question is required. 

In our paper B c -meson decays into heavy quarkonium (J/ip, , 0(2S I ) or Bg) and a set of charged 7r-mesons (wr) ch 
are considered. In the framework of QCD factorization theorem the amplitude of these processes can be splitted 
into two independent parts that describe B c — > V(P)W vertex and W* — > nn transition. The first part is expressed 
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Figure 5: Transferred momentum distributions of B c —± J/ip + (57r) c h (left figure) and B c — > ip(2S) + (57r) c fc (right figure). 
Notations are same as in fig|4] 

(a) (b) 

rf Br 10~ S d Br 10" 5 



dq 2 GeV 2 dq 2 GcV 2 




2 , GeV 2 



Figure 6: Trasferred momentum distributions for B c — >■ By' + (37r) ch (left figure) and B c —> Bl*' + (57r) ch (right figure). Solid 
and dashed curves correspond to fromfactor sets SR and PMf 



through i? c -meson formfactors, that can be calculated using different nonperturbative methods: QCD sum rules, 
various potential models, etc. Information about W* — > 57r transition amplitude, on the other hand, can be obtained 
from theoretical and experimental analysis of r-lepton decays r — > v T + nix. Using this procedure we calculated 
branching fractions and distributions over invariant mass q 2 — M 2 n ^ for n — 3 and 5. It is shown that in the case of 

B c — > J/ip + (37r) c hdecay our predictions agree well with experimental results, presented in (l5| . Other considered in 
our article decays were not observed yet, but their experimental investigation can be expected in the nearest future. 

It should be also noted, that experimental analysis of considered in our article decays could be useful for studying 
hadronization of virtual PF-boson into a set of light mesons. In the ratio 

PT (q 2 ) - (rfBr [B c -s- V(P) + mr] /dq 2 )/{dBr [B c V(P) + iv\ /dq 2 ) (22) 

all dependence on _B c -meson formfactors is canceled, so it describes only W — ¥ mr transition. Previously this function 
was studied only in r-lepton decays, but in _B c -meson decays one can probe higher values of transfered momentum. 
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